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We propose a mechanism of nonthermal production of a neutralino cold dark matter particle, x, from 
the decay of cosmic strings which form from the spontaneous breaking of a U(l) gauge symmetry, 
such as Ub-l(X), in an extension of the minimal supersymmetric standard model (MSSM). By 
explicit calculation, we point out that with a symmetry breaking scale n of around 10 s GeV, the 
decay of cosmic strings can give rise to £l x ~ 1. This gives a new constraint on supersymmetric 
models. For example, the dark matter produced from strings will overdose the universe if rj is 
near the electroweak symmetry breaking scale. To be consistent with Q x < 1, the mass of the 
new (7(1) gauge boson must be much larger than the Fermi scale which makes it unobservable in 
upcoming accelerator experiments. In a supersymmetric model with an extra (7b_l(1) symmetry, 
the requirement of fl x < 1 puts an upper bound on the neutrino mass of about 30eV provided 
neutrino masses are generated by the see-saw mechanism. 



PACS numbers: 98.80Cq 

I. INTRODUCTION 

In spite of the increasing evidence that cold matter 
(matter with pressure p = 0) makes up less than the 
critical density p c for a spatially flat Universe, equally 
strong evidence for the existence of a substantial amount 
of cold dark matter (CDM) remains. The best current 
estimates give flcDM ~ 0.3 |lj whereas fls < 0.1 ||] 
(here, fix = px/Pc denotes the fractional contribution 
of X matter to p Cl and B stands for the contribution of 
baryons) . 

The leading candidates for cold dark matter are the 
axion and the neutralino. The axion is a neutral spin- 
zero Pseudo-Goldstone boson associated with the spon- 
taneous breaking of the global Upq (1) symmetry, which 
was introduced by Peccei and Quinn |3J as a solution to 
the strong CP problem. At zero temperature the axion 
mass is given by 

1f) 12 CleV 

m a ~6x 10~ 6 eVN( W ° eV ) 

J a 

where f a is Peccei-Quinn symmetry breaking scale and 
N is a positive integer which describes the color anomaly 
of Upq(1). Axions can be produced by three different 
mechanisms: vacuum alignment, axion string decay and 
axion domain wall decay Q . Cosmology yields an upper 
limit on f a of /„ < 10 12 GeV. 

The neutralino is an electrically neutral hypothetical 
particle which arises in supersymmetric models. In many 
such models, e.g. in the MSSM (the minimal super- 
symmetric standard model), the lightest supersymmetric 
particle (LSP) is stable, unless R-parity violating inter- 
actions are included. The LSP is generally thought to be 



the lightest neutralino \- The neutralinos in the Universe 
today are in general assumed to be a relic of an initially 
thermal neutralino distribution in the hot early Universe. 
Based on this thermal production mechanism, there have 
been many calculations of the LSP abundance (for a re- 
view, see e.g. ||) as a function of the MSSM parameters. 
These studies show that there exists a domain of parame- 
ter space in the MSSM which is consistent with all of the 
present experimental constraints and for which the x nas 
a relic mass density ft x ~ 1. However, cosmology also im- 
poses limits on the LSP mass. In the case of a Bino-like 
LSP, the calculation of Refs. [§ yields M s < 300 GeV. 
A recent study relaxes this upper bound to about 600 
GeV by including the B coannihilations with the e and 
p. 

In this paper, we propose a new non-thermal produc- 
tion mechanism of the LSP. We consider models with an 
extra U(l) gauge symmetry in extensions of the MSSM. 
This U(l) symmetry could be E/b_l(1), where B and L 
are respectively baryon and lepton numbers. Such mod- 
els explain the neutrino masses via the see-saw mecha- 
nism. Another possibility is that the new U(l) corre- 
sponds to a U(l)' from string theory or grand unified 
theories M. 

The basic idea of our mechanism is as follows. When 
the extra U(l) symmetry which we have introduced gets 
broken at a scale r\, a network of strings is produced by 
the usual Kibble mechanism Q . The initial separation of 
the strings is microscopic, of the order A -1 ?? -1 (where A 
is a typical Higgs self coupling constant of the U(l) sector 
of the theory) which implies that a substantial fraction of 
the energy density of the Universe is trapped in strings. 
After the symmetry breaking phase transition, the defect 
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network coarsens. In the process, string loops decay. If, 
as we assume, the fields excited in the strings couple to 
the neutralino Xj then a non-thermal distribution of \ 
particles will be generated during the process of string 
decay. The total energy density in x particles will depend 
on the scale 77 of U(l) symmetry breaking. The presence 
of our alternative generation mechanism for \ particles 
relaxes the constraints on the mass of the x- Even if 
the usual thermal generation mechanism is too weak to 
generate Q x ~ 1, our new non-thermal mechanism may, 
for appropriate values of 77, be able to lead to f2 x ~ 1. In 
fact, we find that if 77 < 10 8 GeV and M x ~ lOOGeV, then 
our mechanism will lead to il x > 1, unless the couplings 
of the U(l) sector to x are small. Note that there are 
similarities between our non-thermal production and the 
mechanism based on preheating proposed in ]To| . 

To begin with, we consider a general case and calculate 
the relic mass density of the LSP, then we will move on 
to a discussion of some implications. 



II. LSP PRODUCTION VIA STRING DECAY 



The total CDM density is Vic dm = ^therm + ^nonth- 
During the temperature interval between T c and the LSP 
freezeout temperature T x , LSPs released by decaying 
comic string loops will thermalise very quickly with the 
surrounding plasma, and hence will contribute to Qtherm, 
which should not sensitively deviate from the value cal- 
culated by the standard method ^,0. However, below 
the LSP freezeout temperature, since the annihilation of 
the LSP is by definition negligible, each CDM particle re- 
leased by cosmic string decays will contribute to £l non th- 
We obviously must have 
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< 1. 



(i) 



This will lead us to a constraint (a lower bound) on the 
cosmic string forming scale. We now calculate ttnonth- 

We assume that the strings evolve in the friction dom- 
inated regime so that the very small scale structure on 
the strings has not formed yet. The network of strings 
can then be described by a single length scale £(t) F|. In 
the friction dominated period, the length scale has 
been shown to scale as [Oh 



Local cosmic strings form at a phase transitions associ- 
ated with the spontaneous symmetry breaking of a gauge 
group G down to a subgroup H of G if the first homotopy 
group of the vacuum manifold 7Ti(§) is nontrivial. We 
suppose the existence of such a phase transition which is 
induced by the vacuum expectation value (vev) of some 
Higgs field $, (|$|) = 77, and takes place at a tempera- 
ture T c with T c ~ 77. The strings are formed by the Higgs 
field $ and some gauge field A of G whose generator is 
broken by the vev of $. We assume that the generator 
of G associated with A is diagonal so that the strings are 
abelian. The mass per unit length of the strings is given 
by /J = v 2 - 

During the phase transition, a network of strings forms, 
consisting of both infinite strings and cosmic string loops. 
After the transition, the infinite string network coarsens 
and more loops form from the intercommuting of infi- 
nite strings. Cosmic string loops loose their energy by 
emitting gravitational radiation. When the radius of a 
loop becomes of the order of the string width, the loop 
releases its final energy into a burst of <& and A particles. 
Those particles subsequently decay into LSP, which we 
denote by Xi with branching ratios e and e'. For simplic- 
ity we now assume that all the final string energy goes 
into particles. A single decaying cosmic string loop 
thus releases N ~ 27rA _1 e LSPs which we take to have a 
monochromatic distribution with energy E ~ 

In such scenarios, we thus have two sources of cold dark 
matter which will contribute to the matter density of the 
universe. We have CDM which comes from the standard 
scenario of thermal production; it gives a contribution 
to the matter density n t herm- And we also have non- 
thermal production of CDM which comes from the decay 
of cosmic string loops and gives a contribution £l n onth- 



i(/) = i(u(i 
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where £(t c ) ~ (Ar;) - 1 and A is the Higgs self quartic 
coupling constant. The number density of cosmic string 
loops created per unit of time is given by pT| : 
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dt 



dt 



(3) 



where v is a constant of order 1. We are interested in 
loops decaying below T x . 

The number density of LSP released from ti sp till today 
is given by: 



nt™ th {t Q )=Nv jj (J-) < ; ,/< 



(4) 



where the subscript refers to parameters which are 
evaluated today. £p — £(tp) where tp is the time at 
which cosmic string loops which are decaying at time 
t x (associated with the LSP freezeout temperature T x ) 
have formed. Now the loop's average radius shrinks 



*The friction dominated regime lasts from the time t c at 
which the strings network forms until a time t* ~ (Gfi)~ 1 t c , 
where G is Newton's constant Jll]]. In our scenario, the CDM 
is produced at and below the LSP freezeout temperature T x ~ 
10 2 -10 3 GeV. Hence for T c < (10 10 5 - 10 n )GeV = T c *, when 
the temperature of the Universe reaches T x , the strings are 
still in the friction dominated regime. Since we are looking 
for a lower bound T c on the scale 77 of the strings, and since 
as we will show below T c <SC T*, the time interval of interest 
in our scenario is in the friction dominated regime. 
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at a rate @ = -Ti oops Gfi, where r ;oops is a nu- 
merical factor ~ 10 — 20. Since loops form at time 

tp with an average radius R(tp) — X G^M^tp, they 

have shrunk to a point at the time t ~ A rj a M^iJ,. 

Thus t F - (Ar)f oops M^^|- Now the entropy density 
is s = 2 ^-g*T 3 where 5* counts the number of massless 
degrees of freedom in the corresponding phase. The time 

t and temperature T are related by t = 0.3c/* 2 (T)^r- 
where M p / is the Planck mass. Thus using Eqs.(||) and 
([|) , we find that the LSP number density today released 
by decaying cosmic string loops is given by: 

nonth „ -r -a , T 4 

v nonth _ ' L lsp _ °-'° r ,,\2 r -2 — f -yrl ± X 
Y LSP ~ " - ~^~ eVA 1 Ioo P s9*t c 9*t x M pl ^6 ' 

(5) 

where the subscript on g* refers to the time when g* is 
evaluated. 

The LSP relic abundance is related to Y x by: 

n x h 2 « M x y x s(t )pc(to) _1 /i 2 

pa 2.82 x 10 8 F x tot (M x /GeV) (6) 

where h is the Hubble parameter and M x is the LSP 
mass. Now Y[f p = Y* herm + Y™ nth \ hence by setting 
h = 0.70, Eqs. (||) and ([!]) lead to the following con- 
straint: 

5.75 x 10 8 Y™ nth (M x /GeV) < 1. (7) 

We thus see that Eqs. (||) and ^ lead to a lower bound 
on the cosmic string forming temperature T c . 

Recent measurements of cosmological parameters from 
the cosmic microwave background radiation combined 
with Type IA supernovae show evidence for a cosmologi- 
cal constant. In such a scenario, the relic matter density 
satisfies jl3| JIm^ 2 ~ 0.35. 

In Fig. [N, we have plotted the bound on T c as a func- 
tion of (JM X for both Q x h 2 = 1 and Q x h 2 = 0.35. 
We have set g* Tc = 250, g* Tx = 90, T x = M p i = 
1.22 x 10 19 GeV, and the cosmic string parameters v = 1, 
A = 0.5 and T = 10. The region above each curves corre- 
sponds to n x h 2 < 1 (fl x h 2 < 0.35 respectively), and the 
region below to fl x h 2 > 1 (fl x h 2 > 0.35 respectively); 
this region is excluded by observations. We see that 
if there is a cosmological constant, a slightly stronger 
bound on T c is obtained. 

III. IMPLICATIONS FOR PHENOMENOLOGY 

Our results have important implications for super- 
symmetric extensions of the standard model with ex- 
tra {J(l)'s (or grand unified models with an intermedi- 
ate SU(3) C x SU(2) L x U(1) Y x [/(!)' gauge symmetry). 



Most importantly, the requirement fl non th < 1 imposes 
a new constraint on supersymmetric model building and 
rules out many models with a low scale of a new sym- 
metry breaking which produces defects such as cosmic 
strings. 

Consider, for example, the model with an extra 
l^b-l(I) gauge symmetry. In this model, the spectrum of 
the standard model is extended to include right-handed 
neutrinos iV,. The light neutrinos receive masses via the 
see-saw mechanism and the matter-antimatter asymme- 
try of the universe is generated by the out-of-equilibrium 
decay of these right-handed neutrinos. In the latter case, 
leptogenesis can occur by the decay of cosmic strings as- 
sociated with the spontaneous breaking of the Ub-lO-) 
gauge symmetry [ [L4| . In the supersymmetric version of 
this model, the strings will release not only right-handed 
neutrinos Ni, but also their superpartners iVj. The heavy 
neutrinos Ni and their scalar partners JVj can decay into 
various final states including the LSP. The superpotential 
relevant to the decays is 

W = H x tL Vl E c + H 2 eLy v N c , 

where H\, H2,L, E c and N c are the chiral superfields and 
yi , y„ are Yukawa couplings for the lepton and neutrino 
Dirac masses, mi = yiV\,mn — y„V2, with v\ t 2 being the 
vacuum expectation values of the Higgs fields. At tree 
level, the decay rates of Ni into s-lepton plus Higgsino 
and lepton plus Higgs are the same and they are smaller 
than the rate of N% decaying into s-lepton plus Higgs and 
Higssino plus lepton by a factor of 2. If the neutralino 
is higgsino-like, the LSP arise directly from the decays 
of the N and N. If the neutralino is bino- or photino- 
like, subsequent decays of s-lepton into binos or photinos 
plus leptons will produce the LSP. For reasonable values 
of the parameters, we estimate the branching ratio e of 
the heavy particle decay into LSP to be between 0.1 and 
0.5. From Eq. (||) it follows that string decays can easily 
produce the required amount of LSP. However too many 
LSPs will be generated unless the B — L breaking scale, 
Ab-l is higher than about 10 8 GeV . In turn, this will set 
a lower limit on the neutrino masses generated by the see- 
saw mechanism, m„ ~ m 2 D / Kb-l- Inserting numbers 
and taking mo ~ m T ~ 1.8 GeV, one obtains that m v < 
30 eV. 

In models with spontaneous breaking of a Ub-l{X) 
gauge symmetry, upper and lower bounds on the B — L 
breaking scale have already been derived from consid- 
erations of cosmic rays from string decay |l5f| and from 
leptogenesis Jig] , respectively. Our lower bound on the 
B — L breaking scale is independent of leptogenesis. 

Our lower limit on the B — L symmetry breaking scale 
in gauged B — L models and in general models with an 
extra U(l) || pushes the mass of the new gauge boson 
far above the Fermi scale, rendering it impossible to test 
the new physics signals from the extra Z' in accelerators. 

To summarize, we have pointed out a new production 
mechanism for neutralino dark matter which can be ef- 
fective in many models beyond the MSSM, models with 
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FIG. 1. The critical temperature T c 
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branching ratio and the LSP mass for Sl x h = 1 (solid line) 
and Q. x h 2 = 0.35 (dashed line). The region above the curves 
corresponds to Q. x h 2 < 1 (£l x h 2 < 0.35 respectively) and the 
region below correponds to Q x h 2 > 1 (Q, x h 2 > 0.35). The 
latter is excluded by observations. 



extra gauge symmetries which admit topological defects. 
The decay of these defects gives rise to a nonthermal con- 
tribution to the neutralino density. We have focused on 
the nonthermal production of LSPs from string decays. 
Similarly, one could consider LSP production from other 
topological defects. We have calculated the relic LSP 
mass density £l non th as a function of the string scale, the 
freezeout temperature and the mass of the LSP. The LSP 
mass density has two contributions, one from thermal 
production which has been calculated by many authors in 
the literature before, another is the non-thermal produc- 
tion calculated in this paper. Our results indicate that 
if the scale rj of string production is about 10 8 GeV, then 
our nonthermal mechanism can produce the required clo- 
sure density of LSPs. For values of rj smaller than the 
above bound, the model is in conflict with observations 
since the LSPs would overdose the Universe. 

One important caveat must be made concerning our 
calculations. Cosmic strings arising in supersymmetric 
models are generically superconducting |l7]] . In this case, 
the string dynamics may be very different from that of or- 
dinary strings, the dynamics assumed in this paper, and 
thus the corresponding constraints on particle physics 
model building would be quite different. Nevertheless, 
the main point that cosmic string decay in extensions of 
the MSSM can yield a new production mechanism for 
dark matter remains unaffected. 
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